Dynamics and Scaling of 2D Polymers in a Dilute Solution 
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The breakdown of dynamical scaling for a dilute polymer solution in 2D has been suggested by 
Shannon and Choy [Phys. Rev. Lett. 79, 1455 (1997)]. However, we show here both numerically and 
analytically that dynamical scaling holds when the finite-size dependence of the relevant dynamical 
quantities is properly taken into account. We carry out large-scale simulations in 2D for a polymer 
chain in a good solvent with full hydrodynamic interactions to verify dynamical scaling. This is 
achieved by novel mesoscopic simulation techniques. 
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The dynamics of polymer chains has attracted atten- 
tion for decades already. In 3D, polymer dynamics ex- 
hibits rich and complex behavior which depends on the 
solvent conditions and polymer concentration Q . The 
2D case, however, has attracted much less attention. Re- 
cently, it has been realized that it has important applica- 
tions in the field of colloids and biomolecules. Examples 
include the 2D diffusion of DNA oligonucleotides con- 
fined to interfaces [3j and the lateral diffusion of lipids 
and proteins along biological interfaces |4( such as cell 
membranes. Further, the dynamics of polymers in 2D is 
of major importance in thin films whose thickness is less 
than the size of the polymer. Wetting, surface adhesion, 
and flow in confined geometries are examples [j| of this 
broad and fundamental field. 

An important feature of essentially all the 2D diffusion 
processes in soft matter is that they take place in a sol- 
vated environment, which implies that the role of the hy- 
drodynamic interaction (HI) cannot be disregarded. The 
HI originates from interactions mediated by the solvent in 
the presence of momentum conservation. Hydrodynam- 
ics plays a major role in a wide range of applications, and 
hence the understanding of its impact on the dynamics 
of polymer systems has received a great deal of attention. 
In 3D the effects of hydrodynamics are well understood: 
it is well known that the dynamics of polymers in dilute 
solution is well described by the Zimm model 0. In 2D, 
however, the situation becomes significantly more com- 
plicated as will be discussed below. 

To understand the dynamics of polymer chains, with or 
without hydrodynamics, a common technique is to apply 
the theory of dynamical scaling [2j. Two key quantities 
here are the radius of gyration R g and the center-of-mass 
(CM) diffusion coefficient D of the chain. In the dilute 
limit, they follow the scaling relations R g ~ N v and D ~ 
A~" D , with corresponding scaling exponents v and vd, 
respectively. Another central quantity is the intermediate 



scattering function defined as 

S(k,t) = (l/iV)^(exp{ifc- [f m (t) - r n (0)]}) , (1) 

where N is the degree of polymerization, k is a wave 
vector, and {r ra }'s are the positions of the monomers. 
This function should then scale as Q 

S{k,t) = k- 1/u F{tk x ), (2) 

where x is the dynamical scaling exponent related to the 
other exponents through the relation 

x = 2 + v D /v. (3) 

This is valid for k 6 (2ir/R g , 2n/a), where a is the size of 
a monomer. Equations (0) and © are the cornerstones 
of dynamical scaling of polymers. 

In the purely dissipative case, the values of the scaling 
exponents for polymer chains are well understood 0] . In 
the dilute limit the simple Rouse model gives v = 1/2 
and vd = 1. When proper volume exclusion is taken 
into account, v = 3/4 in 2D and approximately 3/5 in 
3D, while vry = 1 still holds for dilute 3D systems and 
for all polymer concentrations in 2D • 

However, when the HI is taken into account, the situa- 
tion becomes dramatically different. While in 3D theory 
and numerical simulations agree with the prediction of 
the Zimm equations that v = Vd [i.e. x = 3) |t1 15 Ifj UtT] . 
in 2D the situation is less clear. What has been es- 
tablished both theoretically Q, Q] and computationally 
[lllll^l is that in good solvent conditions v = 3/ A in 2D, 
as in the case of no HI. 

The situation with vd is more subtle, however. Us- 
ing lattice-gas simulations Vianney et al. [l2T ] found a 
large positive value of vry = 0.78 ± 0.05. The molecular 
dynamics (MD) simulations of S(k,t) by Shannon and 
Choy 11], in turn, gave x = 2 which would imply that 
vd = 0, if Eq. (PI holds. However, from their MD data 
for D vs. N they concluded that ^d > thus contradict- 
ing the scaling law. They also solved the Zimm equations 
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numerically in 2D and verified the result that x = 2, but 
found that now vn < ^lj ■ These results prompted the 
authors of Ref. jljj to suggest that dynamical scaling is 
broken for 2D polymers. Essentially the very basis of 
polymer dynamics is being questioned. 

In this letter, our objective is to determine the validity 
of dynamical scaling for 2D polymers. To this end, we 
first present analytic arguments which show that when 
finite-size effects are properly taken into account, the 
scaling of D with respect to N is truly logarithmic, lead- 
ing to vd = and thus to x = 2. Following this, we 
carefully extract the exponents v, v>d, and x through ex- 
tensive mesoscopic simulations of a 2D polymer in a good 
solvent with the full HI included. Our results verify both 
that x = 2 and the predicted logarithmic scaling of D, 
and thus we conclude that dynamical scaling is obeyed 
in 2D. 

To overcome the significant difficulties in simulating 
polymers with full hydrodynamic interactions, we em- 
ploy a novel mesoscopic simulation method introduced 
recently by Malevanets and Kapral (MK) |i3IIl| . The 
MK method is essentially a hybrid molecular dynam- 
ics scheme, where the polymer chain is treated micro- 
scopically while the solvent obeys coarse-grained dynam- 
ics. In practice this idea is implemented by choosing the 
monomer-monomer and monomer-solvent interactions as 
in MD simulations, while the conservative interactions 
between the solvent particles are absent as in the ideal 
gas. This description preserves the hydrodynamic modes 
through so-called collision rules. Further, it allows for a 
major speedup compared to other simulation techniques 
such as MD. 

To describe the dynamics of the coarse-grained solvent, 
time is partitioned into segments r and the simulation 
box is divided into collision volumes or cells. The effec- 
tive interactions between the solvent molecules take place 
at each r: this is called a collision event. In a collision 
the velocities of the solvent particles are transformed ac- 
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Here Vi is the 



cording to vi (t + t) 

velocity of the particle i, V is the average velocity of all 
the particles in the cell the particle i belongs to, and uj 
is a random rotation matrix chosen for that particular 
cell. It can be shown |13| that this multiparticle collision 
dynamics conserves the momentum and energy in each 
collision volume, and thus gives a correct description of 
the hydrodynamics of the velocity field. 

Our model system consists of a polymer chain with 
N monomers immersed in a 2D coarse-grained solvent. 
The mass of a solvent particle is set to in, and the 
monomer mass is 2m. The monomer-monomer and 
monomer-solvent interactions are described by a trun- 
cated Lennard- Jones (LJ) potential: 



U LJ {r) = {l e[[alr)12 ~ [alrf]+ ' 



r < 2 1 / 6 cr- 

r > 2 1/6 o\ (4) 
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FIG. 1: The scaling of S(k,t) for a polymer with N = 40 in 
a simulation box with L = 120. Here k £ [1.0, 2.4]. 



Here a and e together with m define the LJ unit system, 
where the unit of time is defined as tlj = o\J m/e. In 
addition to the LJ potential, there is an attractive FENE 
potential between the nearest-neighbor monomers: 



U c (r) = -(aRt/2)ln(l-r 2 /Rl) 



(5) 



where a ~ 7 ea~ 2 and Rq = 2 a. 

The solvent density was set to p = 0.581 <7~ 2 , and the 
simulations were carried out at a temperature fc^T = 
1.2 e, yielding good solvent conditions. The equations 
of motion were integrated using the velocity Verlet al- 
gorithm with a time step St = 0.005 t^j. The choice 
of the parameters that determine the collision dynam- 
ics fixes the properties of the coarse-grained solvent, e.g. 
its viscosity. Here we set the collision time to r = tlj 
and the linear size of the collision volume to l c — 2 a. 
The random rotation angles were chosen from a uniform 
distribution in [0, 27r). The size of the polymer chain N 
varies from 20 to 80 monomers, and the linear system 
size L ranges from 40 a up to 420 a. Periodic boundary 
conditions were employed for all system sizes. The CM 
diffusion coefficient D was determined using the memory 
expansion method presented in Ref. jl5l | . 

First, we checked the scaling of R g with N e [20, 80] 
[lij . Our estimate for the scaling exponent of the radius 
of gyration is v — 0.75 ±0.02, in excellent agreement with 
theory. Next, we computed the dynamic structure factor 
S(k, t) which is depicted in Fig. ^ Our data show the 
best collapse with x = 2.0 ±0.1, and we find that x is not 
particularly sensitive to either N or L when N < L (data 
not shown). This confirms the MD results of Shannon 
and Choy |llj , and shows that the 3D Zimm result x = 3 
is indeed invalid in 2D. 

Next, we address the crucial question of the value of the 
scaling exponent vd for the CM diffusion coefficient. In 
the case of the corresponding system in 3D |8| , the finite- 
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size dependence of D is given by D ~ 1/L. Therefore, in 
principle it is easy to determine D for a fixed chain length 
N by running a series of simulations for different values 
of L, and then extrapolating to L — > oo. By repeating 
this procedure for several values of N, the exponent vd 
can be determined. 

However, in the 2D case the finite-size effects are much 
more subtle due to the infinite range of the HI. We have 
carefully calculated D analytically for a 2D polymer in 
a finite system of size L using various approximations, 
including the approaches presented in Refs. pi. Ull \v\. 
In all cases we find 18] that D follows the scaling relation 

D ~ A\nN - Slny, (6) 

where A and B are constants whose values depend on the 
approximations used [lillli^ . This shows that extracting 
i>d in the "traditional" sense in the thermodynamic limit 
L — > oo is no longer possible. 

To numerically study the scaling of D, we deter- 
mined D for each N <E [20, 80] with different values of 
L. For instance, for N = 30 we considered the cases 
L e {60, 90, 120, 150, 180, 210, 240}. For every N, we ex- 
amined the behavior of D as a function of ln(l/L), and 
found that the behavior indeed is linear. Furthermore, 
we can estimate the exponent vjj in terms of effective 
diffusion coefficients in the following way. We chose cut- 
off values L cut e {10 2 , 10 3 , 10 4 , 10 5 , 10 6 }, and extrapo- 
lated a value D(N, L cut ) for each chain size and cutoff. 
If the data complied with Eq. ©, we should, when plot- 
ting D(N,L cu t) vs. In TV, obtain a set of equally spaced 
straight lines. Each line corresponds to a certain cutoff, 
and the lines should all have the same slope A. As can 
be seen in Fig. this indeed holds within the statistical 
uncertainties of our data. 

Most importantly, Fig. confirms the prediction of 
logarithmic scaling of D with N, which means that vd = 
0. To quantify this, we can extract the exponent vd from 
In D(N, L cut ) vs. In AT: for large values of L, we should 
have InD ~ —i/plnN. The results in Fig. show that 
vd decreases steadily with L as it should. For the largest 
L cut studied here, we find vd ~ 0.05 ± 0.05. 

The analysis above reveals the reason for the suggested 
breakdown of scaling in Refs. [ill fl2T]. While the result 
x = 2 is correct, as verified here, the results in the previ- 
ous studies for vd are simply incorrect because the expo- 
nent has been extracted without proper finite-size scaling 
analysis. Thus, we can conclude that dynamical scaling 
holds for 2D polymers with x = 2, v = 3/4, and vd = 0. 

Finally, we wish to discuss the issue of long-time tails 
in 2D diffusion. The velocity autocorrelation function for 
a tagged particle i is defined as 4>{t) = (vi(t + if) ■ Vi(t')} 
and its integral / °° dt</>(i) yields D. In t he p resence of 
the HI it has been shown that <j>(t) ~ i -1 20], or <p(t) ~ 
[i-\/ln(i)] -1 [2l]], which means that rigorously speaking, 
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FIG. 2: The dependence of D on the degree of polymerization 
N for different cutoffs L cut . The lines correspond to L cut £ 
{10 2 , 10 3 , 10 4 , 10 5 , 10 6 } from bottom to top. 



D is not well defined in 2D. This would seem to invalidate 
the present scaling arguments. However, there are several 
ways to resolve this problem. A standard method is to 
view the diffusion coefficients in 2D as time-dependent 
quantities D(t) [2l| . from which one can define effective 
values of D = D(tf) at some finite time tf. Alternatively 
one can normalize D(t)'s with an appropriately chosen 
"bare" diffusion coefficient Do(t) (e.g. that of a monomer 
0). Such normalized values D{t) / D$(t) should converge 
to a finite result at long times. 

However, there are many cases where it is in fact very 
difficult to observe long-time tails in either experiments 
|22| or in simulations. Further complications may arise 
from the fact that confined geometries can strongly in- 
fluence the asymptotic decay of <fr(t) 23]. In the present 
case, the issue of long-time tails is settled by recognizing 
that in dynamical scaling, the absolute values of the dif- 
fusion coefficients are irrelevant: only the behavior of D 
as a function of the chain or system size matters. Hence, 
we can here use such effective values provided that they 
have been determined in a consistent way. To this end, we 
have simply determined Z)'s over a time interval where 
the coefficients have converged within numerical error. 
More precisely, assuming D ~ £ D /tD, where in is the 
distance over which the chain diffuses during the time 
interval <u, the diffusion coefficients have been measured 
at a point where the chain has diffused a scaled distance 
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FIG. 3: The convergence of D(t) vs. time for several systems 
with different values of N and L. The diffusion coefficients 
have been computed using the memory expansion method (see 
Ref. for details). 



io/Rg = 2 — 4. The convergence of D(t) on these time 
scales is demonstrated in Fig. |3 for several systems with 
different values of N and L. 

In conclusion, we have applied the MK algorithm to a 
dilute 2D polymer solution. The MK method itself has 
proved to be an efficient tool for studies of macromolecu- 
lar systems, especially in the dilute limit where the com- 
putational cost is mainly due to the explicit solvent. The 
technique has enabled us — at a moderate computational 
cost — to study system sizes that have not been previ- 
ously amenable to simulations. This approach together 
with proper finite-size scaling analysis has allowed us to 
solve the controversy regarding the dynamical scaling of 
dilute polymer solutions in 2D with full hydrodynam- 
ics. We have found, in contrast to previous arguments 
[fIIIi2| . that the exponent relation x = 2 + Ud/v is valid 
within numerical error. This justifies the scaling hypoth- 
esis, and shows that the anomalous exponent x — 2 found 
in previous studies is due to the logarithmic scaling of D 
as a function of N. 
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